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Abstract.  Univariate  and  bivariate  procedures  for  investigating  the  properties  of 
single  and  joint  intermittent  stochastic  processes  are  presented.  They  allow  the 
characterization  of  all  the  statistics  of  intermittent  variables  using  a set  of  three 
basic  parameters  in  the  multiffactal  framework,  whatever  the  scales  and  the  in- 
tensity. The  multifractal  formalism  is  then  extended  to  more  than  one  variable 
to  investigate  the  degree  of  dependence  among  random  fields  by  examining  the 
nature  of  their  joint  distribution.  These  formalisms,  which  do  not  require  any 
assumption  about  the  spectrum  or  the  probability  distribution  of  the  data  sets 
under  consideration,  are  finally  illustrated,  first,  by  studying  the  properties  of 
turbulent  velocity  fluctuations  recorded  in  a laboratory  grid-generated  turbu- 
lence experiment  and  in  the  ocean  using  a high-frequency  free-fall  profiler.  Sec- 
ond, we  investigate  the  distributions  of  temperature,  salinity  and  in  vivo  fluores- 
cence (a  proxy  of  phytoplankton  biomass)  time  series  recorded  in  the  coastal 
waters  of  the  Southern  Bight  of  the  North  Sea,  as  well  as  their  potential  correla- 
tion. 


Introduction 

Considering  the  current  awareness  of  the  intermittent 
nature  of  both  physical  and  biological  patterns  and  proc- 
esses in  marine  sciences  ( Gibson , 1991;  Pascual  et  al., 
1995;  Seuront  et  al.,  1996a,  1999;  Jou,  1997;  Jimenez, 
2000),  and  the  emergence  of  hot  topics  such  as  those  re- 
lated to  thin-layer  properties  ( Cowles  et  al.,  1998;  Osborn, 
1998),  there  is  a real  need  to  focus  on  the  precise  nature  of 
both  single  intermittent  variables  and  the  degree  of  correla- 
tions between  joint  intermittent  variables,  especially  those 
related  to  physics  and  biology. 

Recent  theoretical  and  empirical  developments,  con- 
ducted both  on  turbulent  energy  dissipation  rates  and  on 
the  associated  passive  scalars  distributions,  have  demon- 
strated the  multiscale  properties  of  intermittent  processes 
(Pascual  et  al.,  1995;  Seuront  et  al,  1996a,  b,  1999, 
2001a,  b;  Lovejoy  et  al.,  2001).  The  large  numbers  of 
modeling  approaches  devoted  to  the  structure  of  turbulent 
flows  have  only  provided  description  at  a very  limited 
range  of  scales.  Indeed,  neither  Gaussian  diffusion  models 
( Visser , 1997)  nor  direct  numerical  simulations  ( Yamazaki 
et  al.,  1991)  possess  a way  to  change  the  scale  upward  or 


downward.  In  contrast,  scale  change  is  natural  for  multi- 
fractal  processes  ( Seuront  et  al.,  1999).  Moreover,  Gaus- 
sian processes  are  determined  by  only  two  moments,  while 
the  first  three  moments  of  empirical  processes  are  usually 
provided  to  characterize  them.  In  the  statistical  framework 
of  multifractals,  we  propose  a more  global  approach  since 
we  provide  all  the  moments  (even  non-integers),  giving  as 
much  statistical  information  as  the  study  of  probability 
density  directly. 

The  procedures  used  to  test  for  independence  between 
two  processes,  based  on  second  order  statistics  analysis 
(i.e.,  correlation),  become  fallacious  for  intermittent  non- 
Gaussian  fields.  It  is  then  necessary  to  test  for  independ- 
ence by  studying  joint  moments  of  all  orders.  In  addition, 
more  recent  procedures  based  on  examination  of  probabil- 
ity density  functions  do  not  deal  with  the  multiscaling 
properties  of  intermittent  fields. 

In  this  paper,  first  we  define  the  concept  of  intermit- 
tency  and  then  review  briefly  the  multifractal  framework. 
Second,  we  will  introduce  an  innovative  objective  tech- 
nique for  determining  if  two  stochastic  processes  can  be 
regarded  as  being  independent  or  not,  and  for  investigating 
the  nature  of  their  potential  coupling.  Finally,  we  provide 
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some  illustrative  case  studies  demonstrating  how  these 
multifractal  formalisms  can  improve  our  understanding  of 
the  structures  and  functions  in  marine  systems. 


£,,(*)~rA 

(1) 

EP(k )«*-*' 

(2) 

Intermittency 

The  concept  of  intermittency  has  been  initially  de- 
scribed in  the  framework  of  turbulence,  where  experimen- 
tal records  of  the  kinetic  energy  dissipation  rate  £ showed 
strong  and  unpredictable  bursts  instead  of  a steady  behav- 
ior of  the  fluctuations  (Fig.  1).  It  has  thus  been  known  for  a 
long  time  that  the  assumption  of  a spatially  uniform  turbu- 
lent cascade  used  by  Kolmogorov  (1941)  is  not  satisfied 
for  the  increments  of  turbulent  velocity  at  distances  in  the 
inertial  subrange  or  shorter,  and  that  strong  gradients  are 
more  common  than  they  would  be  for  a Gaussian  distribu- 
tion ( Batchelor  and  Townsend , 1949).  This  means  that 
most  of  the  time  the  gradients  would  still  be  of  the  order  of 
their  standard  deviation,  but  that  occasionally  we  should 
expect  stronger  burst,  much  more  often  than  in  the  Gaus- 
sian case.  This  is  the  phenomenon  of  intermittency,  which, 
as  with  other  turbulent  phenomena,  is  found  at  all  scales. 

This  may  be  understood  from  a geometrical  point  of 
view  considering  that  the  turbulent  regions  do  not  fill  the 
whole  space,  but  only  a subpart,  in  a very  irregular  way. 
This  erroneous  picture  of  a non-space  filling  turbulence 
takes  the  intermittent  nature  of  turbulence  into  account  by 
assuming  that  subeddies  in  the  inertial  subrange  are  either 
‘dead’  (inactive)  or  ‘alive’  (active),  and  leads  to  a (mono-) 
fractal  description  of  turbulence.  However,  it  is  now 
known  that  both  turbulent  velocity  and  passive  scalars 
fluctuations  are  intermittent  in  the  sense  that  strong  (i.e. 
‘more  active’)  subeddies  occupy  tiny  fractions  of  the  space 
available.  This  led  to  a multifractal  description  of  turbu- 
lence, each  intermittency  level  being  associated  with  its 
own  fractal  dimension  [see  Frisch  (1995)  for  more  details]. 

In  the  following  we  quickly  review  how  to  characterize 
the  statistics  of  a multifractal  intermittent  field. 

Investigating  single  intermittent  distributions: 
a brief  review  of  the  multifractal  framework 

Spectral  analysis  has  been  widely  used  in  ecology  to 
separate  and  measure  the  amount  of  variability  occurring 
at  different  wavenumbers.  In  particular,  assuming  local 
isotropy  and  three-dimensional  homogeneity  of  turbulence 
in  the  inertial  subrange,  the  velocity  fluctuations  and  the 
fluctuations  of  a passive  scalar  quantity  P can  be  de- 
scribed in  Fourier  space  by  the  spectral  densities  Ev  ( k ) 
and  Ep  ( k ) : 


where  k is  a wavenumber,  and  with  /?,,  = /),,  =5/3  for  a 
non-intermittent  turbulence  ( Kolmogorov , 1941;  Obukhov , 
1941,1949;  Corrsin,  1951).  However,  spectral  analysis  and 
related  techniques  are  implicitly  based  on  Gaussian  hy- 
pothesis untenable  in  an  intermittent  framework  ( Seuront 
et  al.,  1999,  2001a).  Thus,  while  random  variability  has 
been  often  modeled  in  marine  ecology  in  the  Gaussian 
framework,  through,  e.g.,  Gaussian  distributions  and 
Brownian  motion  (e.g.,  Yamasaki  and  Okubo,  1995;  Kis- 
ser, 1997),  we  will  here  generalize  this  approach  fully  tak- 
ing into  account  the  intermittency  of  turbulent  velocity  and 
passive  scalar  fluctuations. 

Under  fairly  general  conditions,  the  probability  distribu- 
tion of  a random  variable  X is  equivalently  specified  by 
its  statistical  moments.  Here  we  consider  the  moments  of 
the  increments  of  a random  field  at  all  scales;  this  corre- 
sponds to  the  introduction  of  the  structure  functions  expo- 
nents (q)  and  C,P(q) , which  describe  the  multiscaling 
of  the  statistical  moments  of  order  q of  turbulent  velocity 
and  passive  scalar  fluctuations: 

((AK,  )’)«/<-'*>  (3) 


((AP,y)~i;"ui)  (4) 

where  / is  the  spatial  scale  and  the  angle  brackets  “ Q ” 

indicate  ensemble  averaging.  For  monoscaling  (i.e.,  mon- 
ofractal)  processes,  the  function  £(q)  is  linear: 
C{q)  = q/2  for  Brownian  motion,  and  £(<?)  = <7/ 3 for 
Kolmogorov-Obukhov-1941.  For  multiscaling  processes, 
this  exponent  is  non-linear  and  concave. 

While  the  functions  £v(q)  and  £r(q)  could  depend  on  a 
very  large  number  of  parameters  (a  priori  an  infinite  num- 
ber of  both  integer  and  non-integer  values  of  q ),  in  the 

framework  of  universal  multifractals  they  can  be  deter- 
mined by  only  three  parameters,  with  the  general  form 
C(q)=  Aq  + Bqa  ( Schertzer  and  Lovejoy,  1 987).  For  turbu- 
lent velocity,  the  normalization  conditions  corresponding 
to  the  conservation  of  turbulent  kinetic  energy  within  the 
inertial  subrange  give  (3)  = 1 , so  that  the  resulting  func- 
tion can  be  written  as 


(5) 


For  a passive  scalar  P , there  is  no  such  a relation,  and  we 
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must  keep  three  parameters,  writing  the  expression  as 
(Seuront  et  al.,  1996a,  b,  1999,  2001a) 

£,(*)  = -^(<7"' -<?)•  (6) 

H is  the  degree  of  non-conservation  of  the  average  proc- 
ess: H = 0 for  a conservative  process  (i.e.,  scale- 

independent)  and  H*0  for  a non-conservative  process 
(i.e.,  scale-dependent).  H is  given  by  H =£P(  1),  while  it 
can  be  seen  from  Eq.  (5)  that  £r(l)>l/3  . C,  is  the  codi- 
mension that  characterizes  the  intermittency  of  the  process 
and  satisfies  0 < C,  < d ( d is  the  dimension  of  the  obser- 
vation space):  C,  = 0 for  a homogeneous  process,  and  C, 
increases  with  the  intermittency  of  the  process,  indicating 
that  the  field  values  corresponding  to  any  given  level  of 
variability  are  more  scarce.  The  Levy  index  a is  the  de- 
gree of  multifractality,  bounded  between  a = 0 and 
a = 2 , corresponding  to  the  monofractal  case  and  to  the 
maximum,  or  log-normal,  multifractal  case,  respectively. 

In  other  words,  the  second  terms  of  Eqs.  (5)  and  (6) 
express  the  multifractal  intermittent  deviation  from  mon- 
ofractality,  in  which  case  £v(,q)  - q/3  and  £P(q)  = qH  . 
The  “distance”  between  monoffactality  and  multifractality 
is  then  a function  of  C,  and  a . The  knowledge  of  these 
parameters  is  enough  to  characterize  all  the  one-point  sta- 
tistics of  the  fields.  In  particular,  we  see  from  Eqs.  (5)  and 
(6)  that  C,  is  given  by  Cl£  = l-3£j(3)  for  velocity  fluctua- 
tions, and  Cu  - H -£P(l)  for  a passive  scalar  quantity  P . 
a and  a are  estimated  as  the  best  nonlinear  fit  of  Eqs.  (5) 
and  (6)  for  values  between  0 and  2,  using  a simplex  proce- 
dure. 

Finally,  we  note  here  that  in  the  multiscaling  frame- 
work, intermittency  is  taken  into  account  noting  that 

& =1  + ^(2)  (7) 

Pr  =l  + C/>(2)  • (8) 

We  thus  see  that  the  intermittency  corrections  introduced 
by  the  second  term  of  Eq.  (5)  lead  to  pr  > 5/3  . Further- 
more, in  order  to  conduct  direct  quantitative  comparisons 
between  the  velocity  and  passive  scalar  intermittent  correc- 
tions associated  to  C,  values,  let  us  note  that  the  coeffi- 
cients in  front  of  the  non-linear  terms  in  Eqs.  (5)  and  (6) 
are  similar  only  if  we  compare  Clv  and  Clp , with  Cxv 
estimated  as  ( Schmitt  et  al,  1996;  Seuront  et  al.,  1996a) 

C\v  - 3"°'  Cte . (9) 


Investigating  joint  intermittent  distributions: 
the  “Generalized  Correlation  Functions” 

Previous  standard  procedures  devoted  to  test  for  inde- 
pendence between  two  given  processes  were  generally 
based  on  second  order  statistics  (i.e.,  covariance  and  corre- 
lation functions),  even  when  they  were  conducted  in  scal- 
ing framework  related  to  spectral  analysis  ( Legendre  and 
Legendre,  1984)  or  geostatistical  analysis  ( Kitanidis , 
1997).  More  recent  procedures  are  based  on  probability 
density  functions  examination  ( Lueck  and  Wolk,  1999). 
The  former  are  implicitly  based  on  Gaussian  hypothesis 
(uncorrelation  implying  independence),  untenable  in  an 
intermittent  framework  characterized  by  a high-order  sta- 
tistical behavior.  The  latter  do  not  deal  with  the  intrinsic 
multiscaling  properties  of  intermittent  fields.  We  then  pro- 
pose here  an  original  testing  procedure  based  on  a high- 
order  generalization  of  the  correlation  concept  between 
two  variables  X and  Y . 

Instead  of  random  variables  X and  Y , we  consider 
here  two  stochastic  processes  AX,  and  AY,  (e.g.,  Parzen, 

1962).  For  more  convenience,  let  us  note  x = AX,  = e°'  and 
y = AY,  = ec- . The  joint  moments  can  be  written  as  the 
moments  of  a vectorial  process: 

(//)  = (e"°'+’C;  ^ ^ x /5(d)  (10) 

where  the  vectors  Q and  G are  given  by  Q = {p,q)  and 
G = (G,,G2),  and  the  exponents  S(Q)  characterize  the 
multiscaling  properties  of  the  joint  moments  (*V)  • The 

“Generalized  Correlation  Function”  (GCF  hereafter)  we 
introduce  is  then  simply  a normalization  of  the  joint  mo- 
ments following: 

The  “Generalized  Correlation  Exponent”  (GCE  hereafter), 
estimated  as  the  slope  of  the  linear  trend  of  c(p,  q)  vs.  / 
in  a log-log  plot,  is  then  expressed  as 

r(p,q)  = CAP)+U<l)-S(p,q)  (12) 

where  Cx(p)  and  L,y{q)  characterize  the  multiscaling 
properties  of  the  single  fluctuations  ( x pj  and  as  de- 
fined in  Eqs.  (3)  and  (4),  and  S(p,q)  characterizes  the 
multiscaling  properties  of  the  joint  fluctuations  (x^y^ , 
see  Eq.  (11).  Both  c(p,q)  and  r(p,q)  are  generalization  of 
correlation  functions.  They  express  the  correlation  between 
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x’’  and  yq  and  their  scale  and  moment  dependence.  In  the 
particular  case  p = q- 1,  Eq.  (11)  recovers  the  classical 
expression  of  the  correlation  coefficient  between  x and 
y . We  nevertheless  need  to  emphasize  here  that,  whereas 
independence  implies  uncorrelation,  uncorrelation  does  not 
imply  independence. 

Indeed,  uncorrelation  corresponds  simply  to  the  relation 
r(l,l)  = 0.  Uncorrelation  implies  independence  only  in 
special  cases  such  as  for  Gaussian  processes.  To  show  this, 
let  us  consider  the  joint  scaling  function  for  lognormal 
multifractals  x and  y.  Using  results  for  multivariate 

Gaussian  processes  (see  any  text  book  on  multivariate  sto- 
chastic processes,  e.g.,  Samorodnitsky  and  Taqqu,  1994), 
one  has  the  general  expression  for  a lognormal  process: 

S(P,q)  = alp  + a2q-aiq2 -aAp2-(7pq  (13) 

so  that 

r(p,q)  = S(p,0)  + S(0,q)-S(p,q)  (14) 

giving 

r(p,q)=Opq.  (15) 

In  this  case,  it  is  clear  that  r(l,l)  or  r( 2,2)  is  enough  to 

estimate  the  only  parameter,  namely  the  correlation  coeffi- 
cient a,  so  that  if  r(l,l)  = 0 or  r(2,2)  = 0 it  can  be  con- 
cluded that  the  two  processes  are  independent.  In  the  gen- 
eral case,  this  is  no  more  true:  independence  between  the 
stochastic  processes  x and  y means  that  the  GCE  verifies 
r(p,q)  = 0 whatever  the  values  of  p and  q , while  uncor- 
relation corresponds  to  r(l,l)  = 0. 

A function  related  to  our  c(p,q)  has  also  been  inde- 
pendently introduced  in  the  lognormal  multifractal  frame- 
work for  financial  time  series  by  Muzy  et  al.  (2000),  where 
it  was  called  “joint  scaling  exponent”. 

The  function  c{p,q ) and  its  related  scaling  exponent 
r{p,q)  can  be  used  as  an  analysis  tool  to  study  some  com- 
plex relationship  between  two  multifractal  fields  x and  y . 
To  provide  some  basis  for  discussion  and  interpretation  of 
experimental  results,  let  us  consider  some  limit  cases. 

If  x and  y are  independent,  as  was  said  above, 
r(p,q)  = 0.  On  the  other  hand,  in  case  of  perfect  propor- 
tionality x = K.y , where  K is  a constant,  or  for  “random 
proportionality”  x — K._y , where  K is  a random  variable 
independent  on  y , it  is  readily  seen  that 

r(p,q)=Cy(p)+Cr(<i)-Cr(p+q)-  (16) 

In  particular,  one  may  note  that  r(p,q)  > 0 due  to  the  con- 
vexity of  the  scaling  functions  C(p) . This  relation  can  be 


directly  tested  to  verify  the  random  proportionality  hy- 
pothesis. Furthermore,  the  shape  of  the  surface  obtained  is 
symmetric  in  the  p-q  plane. 

Other  simple  situations  may  be  considered:  if  x = K.y" 
with  b > 0 and  K constant,  or  if  x-  K.yh  with  K random 
and  independent  of  y , then  one  has 

r(p^)=CAbP)+Cr(q)-Cr(bP+q)  ■ 0?) 

Eq.  (17)  is  still  positive,  but  no  more  symmetric  in  the  p-q 
plane;  it  is  symmetric  in  the  bp-q  plane.  In  this  framework, 
the  value  of  b may  be  first  estimated  as  the  positive  value 
such  that 

r(p,0)  = r(0,bp)  (18) 

Using  the  value  estimated  this  way,  this  framework  is  then 
tested  by  verifying  that  r(p/b,q)  is  indeed  symmetric  in 
the  p-q  plane.  More  generally  speaking,  the  more  r(p,q)  is 
positive,  the  more  the  x = AX,  and  y = AY,  are  dependent 
random  variables. 

The  main  advantages  of  this  framework  are  the  follow- 
ing: (i)  it  makes  no  assumptions  about  the  spectrum  or  the 
distribution  of  either  data  sets;  (ii)  it  takes  fully  into  ac- 
count their  intrinsic  multiscaling  properties;  (iii)  it  may  be 
used  to  provide  much  more  fundamental  relations  between 
two  pattems/processes;  and  iv)  more  generally  it  fully  ex- 
plores qualitatively  and  quantitatively  the  correlations  of 
the  fields  x''  and  yq , revealing  the  correlations  of  large 
and  low  fluctuations  of  both  processes. 

Turbulent  velocity  fields 
Turbulence  measurements 

Laboratory  experiments 

T urbulence  has  been  generated  by  means  of  fixed  PVC 
grids  (diameter  2 mm,  mesh  size  1 cm)  in  a circular  flume. 
Considering  the  dimensions  of  the  flume  and  the  flow  ve- 
locities considered  in  the  experiment  (Table  1),  we  were 
dealing  with  a fully  developed  turbulence.  A high- 
frequency  (100  Hz)  time  series  of  instantaneous  horizontal 
turbulent  velocity  was  measured  by  hot-film  velocimetry 
10  and  20  cm  behind  the  grid.  The  average  turbulent  en- 
ergy dissipation  rate  was  derived  following  the  isotropic 
formula  ( Tennekes  and  Lumley,  1972): 

£ = 1 5v  J*”“  k2E(k)dk  (19) 

where  £ is  the  turbulent  energy  dissipation  rate  (mV3), 
k the  wavenumber  (k  = 2n!X,  m'1),  X the  eddy  wave- 
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length  (m),  kmm  and  k^,  respectively  the  smallest  and 
largest  wavenumber  resolved,  and  E(k)  the  turbulence 
spectrum  (m  s'2).  The  spectrum  E(k)  is  regarded  as  the 
mean-square  amplitude  of  velocity  fluctuations  associated 
with  a wavenumber  of  turbulent  motion. 

Table  1.  Main  characteristics  of  the  grid-generated  turbulence 
experiment  (d:  distance  from  the  grid;  Re:  Reynolds  number;  N: 
numbers  of  time  series  analyzed). 


d (cm) 

Speed  (m  s'1) 

e (m2  s'3) 

Re 

N 

20  cm 

0.15 

2.86  10's 

45  10! 

5 

0.20 

3.54  10-6 

60  103 

5 

0.25 

9.14  10-6 

75  103 

5 

0.30 

1.36  10's 

105  103 

5 

0.35 

1.64  Iff5 

105  103 

5 

0.40 

1.82  10'5 

120  103 

5 

10  cm 

0.15 

1.01  10'7 

45  103 

4 

0.20 

1.24  10'5 

60  103 

5 

0.25 

3.2  10'5 

75  103 

5 

0.30 

4.75  10'5 

105  103 

5 

0.35 

5.75  10'5 

105  103 

5 

0.40 

6.36  10'5 

120  103 

5 

To  investigate  potential  differences  in  the  structure  of 
velocity  fields  in  relation  to  the  different  levels  of  the  tur- 
bulent kinetic  energy,  we  analyzed  here  5 time  series  of  1 
min  duration  (i.e.,  6000  data  points)  for  5 flow  velocities  at 
the  two  locations  considered  behind  the  grid. 

Field  experiments 

Measurements  of  the  turbulence  fields  were  done  using 
TurboMAP  (Turbulence  Ocean  Microstructure  Acquisition 
Profiler),  a free-fall  integrated  system  housed  in  an  HCR 
treated  aluminum  pressure  case,  specifically  designed  for 
concurrent  measurements  of  turbulent  parameters  (shear, 
duldz , where  u and  z are  the  turbulent  fluctuating  cross- 
stream velocity  component  and  the  depth,  respectively), 
hydrographic  parameters  (pressure,  temperature  and  con- 
ductivity) and  bio-optical  parameters  (in  vivo  fluorescence 
and  backscatter).  All  sensors  are  set  into  a parabolically 
shaped  cap  at  the  top  of  the  pressure  case.  This  arrange- 
ment ensures  that  these  sensors  point  into  an  undisturbed 
flow.  All  signals  are  simultaneously  sampled  at  256  Hz 
and  stored  internally  in  a PCMCIA  memory  card.  More 
details  on  the  description  and  evaluation  of  both  design 
and  performance  of  this  instrument  can  be  found  in  Wolk 


et  al.  (2001  a,b). 

In  this  paper,  we  focus  only  on  the  ability  of  this  new 
high-resolution  biophysical  profiler  to  measure  the  shear 
signal,  duldz . In  that  way,  we  use  8 vertical  shear  profiles, 
recorded  on  21  August  1998  in  a tidal  channel  located  be- 
tween Neko-Seto  and  Meneko-Seto  Straits  (Japan).  The 
sampling  location  has  been  chosen  because  this  area  is 
characterized  by  very  strong  tidal  currents  (up  to  180  cm  s" 
'),  associated  with  channel  narrowing,  and  a subsequent 
high  tidal  mixing  ( Takasugi  et  al.,  1994a,  b). 

The  turbulent  velocity  fluctuations  are  measured  with  a 
standard  piezo-ceramic  shear  probe,  which  has  extensively 
been  described  elsewhere  (e.g.,  Lueck  et  al.,  1997).  The 
probe  is  associated  with  three  orthogonal  accelerometers 
mounted  inside  the  pressure  case  and  is  used  to  quantify 
noise  level  of  velocity  shear  and  to  clean  the  shear  signal 
of  possible  narrow-band  contaminations  from  vehicle  vi- 
brations (Levine  and  Lueck,  1999). 

To  recover  the  turbulent  fluctuating  velocity  u from 
duldz,  we  run  an  anti-derivative  digital  filter  over  the 
shear  data,  which  consists  of  a scaled  single-pole  low-pass 
filter  with  a cut-off  frequency  fc . The  fc  is  a function  of 

the  integration  time  scale  tt  as  f(  =1/2.3/,  [see  Wolk  and 
Lueck  (2001)  for  further  details].  Here  we  used  f,  = 7.5  s 
and  f.  = 0.06  Hz,  leading  to  a maximum  wavelength  re- 
solved to  10  m.  The  lowest  frequency  resolved  is  deter- 
mined by  the  high-pass  filtering  of  the  shear  probe  signal 
to  minimize  pyro-electric  effects  (70  Hz). 

The  vertical  distribution  of  the  turbulent  kinetic  energy 
dissipation  rates  £ has  been  estimated,  after  systematic 
denoising  and  despiking,  from  shear  power  spectra  com- 
puted for  512  successive  segments  of  shear  (i.e.,  every 
meter  considering  an  average  profiler  sinking  velocity  of 
0.5  m s'1)  following  the  isotropic  formula  defined  above  in 
Eq.  (19). 

Considering  the  extreme  fluctuations  observed  in  the 
vertical  distribution  of  the  dissipation  rates  £ (Fig.  1),  our 
analysis  has  been  conducted  on  all  the  512  data  point  seg- 
ments available  within  each  vertical  profile  (Table  2). 
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Figure  1.  Example  of  the  intermittent  character  of  the  vertical 
distribution  of  the  kinetic  energy  dissipation  rate  £ , which  has 
been  estimated  every  meter  using  Eq.  (19). 


Table  2.  Depth,  numbers  of  512  data  points  subsections,  and 
mean  of  the  turbulent  kinetic  energy  dissipation  rate  £ (m‘  s'3), 
for  each  of  the  eight  vertical  profiles. 


Data  set 

N 

n 

£ 

1 

14463 

28 

1.80  10* 

2 

20479 

40 

1.87  10'7 

3 

21887 

43 

6.98  10‘8 

4 

18431 

36 

6.43  10'8 

5 

10495 

20 

9.61  10'8 

6 

22911 

45 

1.33  10-6 

7 

29949 

58 

3.72  10'7 

8 

17279 

34 

5.65  10'8 

Characterizing  turbulent  velocity  fields 

The  velocity  spectra  exhibit  a scaling  behavior  ranging 
between  0.1 85 ±0.1  m and  0.013 ±0.015  m for  velocity 
fluctuations  measured  in  our  circular  flume,  and  between 
1.0  ±0.1  m and  0.050  ±0.005  m for  the  velocity  fluctua- 
tions estimated  from  shear  data  recorded  in  the  field  (Fig. 
2).  These  scaling  behaviors  are  qualitatively  close  (we  try 
a quantification  later)  to  the  theoretical  -5/3  Kolmogorov’s 
power  law  over  roughly  1.1  and  1.3  decades,  respectively, 
with  Py  ~ 1 .70  in  both  cases.  The  ranges  of  scales  defined 
that  way  will  be  regarded  as  the  inertial  subrange. 

The  estimates  of  the  velocity  structure  functions 
^( AF confirm  the  scaling  regimes  previously  shown  by 
spectral  analysis  for  different  orders  of  moment  q (Fig.  3). 
Estimated  scale  ratios,  highest  and  lowest  wavenumbers  of 
scaling  ranges  observed  for  the  5 1 2 data  points  subsections 
considered  within  each  profile,  exhibit  very  similar  pat- 
terns to  those  observed  in  the  spectral  framework.  We 


refer  the  reader  to  Table  1 for  further  details.  We  subse- 
quently computed  the  scaling  exponents  £{q)  for  several 
values  of  q between  0 and  10  (i.e.  with  0.1  increments) 
for  the  range  of  scales  where  ^(3)  = 1 is  valid;  this  range 
of  scales  characterizes  here  the  inertial  subrange.  The  scal- 
ing of  the  second-order  moments  £V(2)  confirms  the  esti- 
mates of  the  spectral  exponents  py  [ fiv  = l + f,(2),  Eq. 
(7)]:  £,.(2)  = 0.710±0.010  for  grid-generated  turbulence 
and  £).(2)  = 0.705  ±0.005  for  ocean  turbulence. 


Figure  2.  Samples  of  the  power  spectral  density  E(k)  of 
turbulent  velocity  fluctuations  recorded  with  a hot-film 
velocimeter  (A)  located  behind  a grid  in  a circular  flume, 
and  (B)  in  a tidally  mixed  tidal  channel,  shown  in  log-log 
plots  as  a function  of  the  wavenumber  k (m1).  The  vertical 
dotted  lines  indicate  the  inertial  subrange  over  which  the 
spectra  roughly  follow  the  power-law  form  £(£)<*  , 

with  fa.  =5/3  (continuous  lines). 
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Figure  3.  The  velocity  structure  functions  ([SVl)q'j  of  ocean 

velocity  versus  the  length  scale  / in  a log-log  plot  for  various 
values  of  the  statistical  order  of  moment  q (i.e.,  <7=1,2  and  3, 
from  top  to  bottom).  Linear  trends  are  clearly  visible,  whatever 
the  order  of  moment  q (continuous  lines),  and  indicate  the  mul- 
tiscaling nature  of  our  datasets.  The  vertical  dotted  lines  corre- 
spond to  the  upper  and  lower  limits  of  the  inertial  subrange.  Note 
that  the  different  scaling  ranges  (estimated  from  the  profile 
shown  in  Fig.  2a)  are  indistinguishable  from  the  one  estimated 
from  spectral  analysis. 

More  generally,  the  clear  nonlinearity  of  the  functions 
£v(q)  (Fig.  4)  is  direct  evidence  for  the  multiffactal  nature 
of  the  velocity  fluctuations.  Moreover,  within  experimental 
error,  the  functions  C v(q)  for  grid-generated  and  ocean 
turbulent  velocity  fluctuations  cannot  be  qualitatively 
shown  as  being  different.  While  the  similarity  observed  in 
the  shape  of  the  functions  £y(q)  qualitatively  suggests 
some  convergent  behavior  in  the  high-order  structure  of 
velocity  fluctuations,  the  values  of  the  universal  multifrac- 
tal parameters  C|E  and  a€  lead  to  specifying  these  re- 
sults. These  parameters  have  been  estimated  as 
Cl£  =0.18  ±0.02  , and  ac=  1.52  ±0.03  for  grid-generated 
turbulence,  and  Clc  =0.1 6 ±0.02  and  ac=  1.55 ±0.04  for 
ocean  turbulence.  That  leads  to  confirmation  that  the  struc- 
ture of  velocity  fluctuations  recorded  in  the  field  and  in  the 
laboratory  cannot  be  regarded  as  being  significantly  differ- 
ent. 

On  the  other  hand,  the  parameters  Cw , estimated  using 
Cic  and  ac  in  Eq.  (9)  as  Cxv  =0.032  ±0.02 , appear  slightly 

smaller  than  the  values  reported  in  the  literature  for  the 
parameters  C„,  characterizing  the  intermittent  behavior  of 
temperature  and  salinity  in  highly  dissipative  areas: 


q 


Figure  4.  The  scaling  exponents  structure  function  ^Y(q)  em- 
pirical curves  (open  dots)  estimated  from  grid-generated  turbu- 
lent velocity  fluctuations,  compared  to  the  monofractal  curve 
^ (9)  = <? /3  (dashed  line),  and  to  the  universal  multifractal 
curve  (thick  continuous  curve)  obtained  with  Cu  = 0. 1 8 and 
ae  = 1 .52  inEq.  (5). 

CtTe  [0.035-0.060]  and  Cl5  e [0.035-0.062]  {Seuront  et  al., 

1996a,b,  1999,  2001a).  This  implies  that  the  temperature 
and  salinity  fields  are  much  more  intermittent  than  the 
velocity  field;  this  high  intermittency  can  be  related  to  the 
strong  discontinuities  (“ramps”)  observed  in  the  tempera- 
ture signals  ( Sreenivasan , 1991;  Pumir,  1994b).  Intermit- 
tency  for  passive  scalars,  in  both  two  and  three  dimen- 
sions, has  also  been  studied  in  direct  numerical  simulations 
( Holzer  and  Siggia,  1994;  Pumir,  1994a,b).  Such  simula- 
tions often  reveal  extended  regions  of  almost  uniform  tem- 
perature separated  by  boundary  layers  (ramps)  with  strong 
gradients. 

We  should  finally  note  that  the  multifractal  parameters 
Clc  and  ae  obtained  here  for  the  very  first  time  for  oce- 
anic turbulence  are  very  similar  to  the  parameters  esti- 
mated for  atmospheric  turbulence:  ClE  =0.15  ±0.03  and 
ac=  1.50  ±0.05  (Fig.  5;  Schertzer  et  al.,  1995;  Schmitt  et 
al.,  1992a,  b,  1993,  1996;  Chigirinskaya  et  al.,  1994).  As 
previously  noticed  by  Seuront  et  al.  (2001b),  these  results 
suggest  a strong  similarity  of  the  intermittent  nature  of 
turbulent  velocity  fluctuations  in  the  atmosphere  and  the 
ocean.  The  resolution  of  this  particular  issue,  well  beyond 
the  scope  of  the  present  paper,  is  in  progress.  In  particular, 
the  question  is  whether  atmospheric  and  oceanic  tempera- 
ture fields  behave  in  a similar  fashion  ( Seuront  and  Yama- 
zaki,  200 1 ; Seuront  et  al. , 200 1 c). 
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Figure  5.  Comparison  of  the  empirical  values  of  the  functions 
C,v(q)  estimated  from  grid  generated  turbulence  (open  dots)  and 
field  ocean  turbulence  (open  triangles),  with  the  universal  multi- 
fractal fit  obtained  with  Clt  and  OC£  values  characteristic  of  at- 
mospheric turbulence,  i.e.  Clt=0.15  and  at  =1.50  in  Eq.  (5). 
The  linear  curve  t^l.(q)~ql 3 expected  in  case  of  non- 
intermittent  turbulence  is  shown  for  comparison  (dashed  line). 

Ocean  passive  and  “biologically  active” 
scalars 

Phytoplankton  of  the  Eastern  English  Channel 

In  a recent  paper,  Seuront  et  al.  (1999)  showed  that  the 
small-scale  intermittent  distribution  of  temperature  (a  sca- 
lar quantity  passively  advected  by  turbulent  fluid  motions) 
and  phytoplankton  biomass  recorded  in  the  coastal  waters 
of  the  Eastern  English  Channel  were  indistinguishable, 
both  in  spectral  and  multifractal  framework  (Fig.  6).  They 
subsequently  concluded  their  study  stating  that  phyto- 
plankton biomass  could  be  regarded  as  a purely  passive 
scalar  at  small  scales.  Here,  small  scales  correspond  to 
temporal  scales  bounded  between  1 and  20  seconds.  Using 
the  Taylor’s  theory  of  frozen  turbulence,  these  temporal 
scales  can  be  related  to  spatial  scales  bounded  between  0. 1 
and  22  meters.  For  larger  scales,  the  phytoplankton  bio- 
mass exhibits  a very  specific  heterogeneous  distribution, 
dominated  by  biological  processes  resulting  from  complex 
interactions  between  the  turbulence  level  of  fluid  motions 
(mainly  different  tidal  conditions  related  to  both  ebb/flood 
and  neap/spring  tidal  cycles),  the  phytoplankton  biomass 
concentration,  and  the  specific  composition  of  phytoplank- 
ton assemblages  ( Seuront  et  al.,  1996b;  Seuront  et  al., 
1999).  These  larger  scales  are  not  investigated  in  the  pre- 
sent work. 


Table  3.  Mean  values  of  the  empirical  estimates  of  the  spectral 
exponent  fir  , and  the  universal  multifractal  parameters  H F , 

CIF  and  aF  for  subsets  of  in  vivo  fluorescence  time  series  re- 
corded under  different  levels  of  turbulence. 


£ (m2.s'3) 

P 

H 

c, 

a 

1 to-6 

1.76 

0.43 

0.060 

1.80 

5 lCV6 

1.76 

0.43 

0.055 

1.82 

1 10'5 

1.77 

0.42 

0.050 

1.80 

5 10'5 

1.78 

0.42 

0.045 

1.81 

1 10J 

1.78 

0.41 

0.040 

1.80 

Turbulence  and  phytoplankton  patchiness 

Considering  the  role  played  by  turbulence  in  the  forma- 
tion and  disruption  of  particle  aggregates  [e.g.,  Kiorboe 
(1997)  and  references  therein],  and  the  highly  varying  dis- 
sipation rates  of  turbulent  energy  observed  over  a tidal 
cycle  in  this  area  (between  10'6  and  1 O'4  m2  s'3;  Seuront  et 
al.,  2001a),  the  question  here  is  whether  the  local  hydro- 
dynamic  conditions  have  an  effect  on  both  temperature  and 
phytoplankton  biomass  local  distributions. 

The  original  46  h 24  min  time  series  recorded  at  1 Hz 
(i.e.,  167040  data;  Seuront  et  al.,  1999)  has  then  been  di- 
vided into  46  subseries  of  1 h (i.e.,  3600  data).  We  then 
built  5 sets  of  7 subseries,  each  set  corresponding  to  a dif- 
ferent level  of  dissipation  rates  ranging  from  e = 10'6  m2.s'3 
to  E = 1 O’4  m2.s'3.  We  then  conducted  both  spectral  and 
multiffactal  analyses  on  each  subseries  of  temperature  and 
in  vivo  fluorescence  within  a given  subset,  and  the  result- 
ing observations  led  to  different  comments  (Table  3). 
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Figure  6.  Comparison  of  the  functions  £T(q)  (open  dots) 
and  £F(q)  (black  dots)  estimated  from  small  scale  (<20 
seconds)  temperature  and  in  vivo  fluorescence  fluctuations; 
(adapted  from  Seuront  et  al,  1999). 
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The  slopes  j3r  and  fiF , estimated  for  the  frequency 
bands  that  maximized  the  coefficient  of  determination  and 
minimized  the  total  sum  of  squared-residuals  for  the  re- 
gression in  log-log  plots  of  the  temperature  and  fluores- 
cence power  spectra  [cf.  Eq.  (2)],  cannot  be  statistically 
distinguished  within  each  set  of  temperature  and  fluores- 
cence subseries  (analyses  of  covariance,  p>0.05),  nor 
between  the  5 sets  considered  (p>0.05).  These  results, 
fully  congruent  with  previous  studies  conducted  in  the 
same  area  ( Seuront  et  al.,  1996a,  b,  1999,  2001a),  suggest 
extreme  similarity  between  the  single  scaling  properties  of 
small-scale  temperature  and  fluorescence  distributions.  At 
this  stage,  phytoplankton  biomass  can  then  still  be  re- 
garded as  being  a passive  scalar. 

The  nonlinear  behavior  of  the  functions  £(q)  confirms 
and  generalizes  in  different  hydrodynamic  regimes  the 
multifractal  nature  of  temperature  and  in  vivo  fluorescence 
observed  by  Seuront  et  al.  (1999)  using  a global  approach 
of  the  same  data  set  (Fig.  6).  However,  the  functions  £r(q) 
obtained  for  temperature  time  series  remain  similar  what- 
ever the  external  hydrodynamic  forcing  (i.e.,  the  values  of 
the  dissipation  rate  £ ),  while  the  nonlinearity  (i.e.,  con- 
vexity) of  the  functions  £F(q)  characterizing  phytoplank- 
ton biomass  clearly  increases  when  the  value  of  the  dissi- 
pation rate  £ decreases  (Fig.  7). 
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Figure  7.  Comparison  of  the  empirical  values  of  the  functions 
C,(q)  estimated  for  temperature  (continuous  curve)  and  in  vivo 
fluorescence  for  different  values  of  the  turbulent  kinetic  energy 
dissipation  rate  £ (£  = 1(T6  m2  s'3,  rhombs;  £ = 5xlCT<  m2  s'3, 
circles;  £ = 10'  m2s'\  squares;  £ = 5 x 1 0”5  m2  s"3,  crosses; 
£ = 10^  nr  s'3,  triangles).  The  linear  curve  (r(q)  = qH  ex- 
pected in  case  of  non-intermittent  turbulence  in  shown  for  com- 
parison (dashed  line). 


The  distribution  of  temperature  then  remains  the  same 
whatever  the  hydrodynamic  conditions.  This  result  is  con- 
vergent with  our  previous  observations  regarding  the  dis- 
tribution of  turbulent  velocity  fields,  and  it  confirms  the 
passive  scalar  behavior  of  temperature. 

On  the  other  hand,  the  distribution  of  phytoplankton 
biomass  is  less  intermittent  and  is  closer  to  that  of  tem- 
perature in  high  turbulent  conditions  (Fig.  7).  This  sug- 
gests that  phytoplankton  cells  behave  more  like  a passive 
scalar  under  high  hydrodynamic  conditions  but  exhibit 
very  specific  properties  under  lower  levels  of  turbulence. 

The  universal  multifractal  parameters  H , C, , and  a 

confirm  these  qualitative  observations  (Table  3).  Let  us 
remember  here  that  the  parameters  H , C, , and  a have 
been  obtained  by  fitting  the  universal  multifractal  model 
detailed  in  Eq.  (6)  to  the  empirical  function  'Q(q)  as  de- 
scribed above.  These  parameters  remain  constant  for  tem- 
perature time  series  with  Hr  = 0.40 , Cir  = 0.05  and 

aT  = 1 .90 . On  the  other  hand,  it  is  shown  that  only  H and 
C,  remain  roughly  constant  for  fluorescence,  with 
Hf  =0.42 ±0.01  and  aF  = ] .8 1 ± 0.01  [values  compatible 
with  the  HF  = 0.43  and  aF  = 1 .80  found  by  Seuront  et  al. 
(1999)].  However,  the  parameter  C{F  is  highly  variable.  It 
increases  33%  when  the  dissipation  rate  £ decreases  from 
the  highest  to  the  lowest  values  experienced  over  a tidal 
cycle,  i.e.  from  £ = 1 O'4  to  £ = 1CT6  m2.s°.  Under  low  hy- 
drodynamic conditions  the  higher  values  of  Cif  indicate 
the  occurrence  of  few  patches  of  high  phytoplankton  con- 
centrations that  are  several  orders  of  magnitude  above 
background  levels.  Under  higher  hydrodynamic  conditions 
the  lower  values  of  CiF  indicate  that  phytoplankton  cells 
are  more  homogeneously  distributed  and  behave  as  a pas- 
sive scalar.  Here  the  so-called  homogenization  effect  is 
thought  to  be  associated  with  the  disruption  of  these 
phytoplankton  patches.  This  hypothesis  is  supported  by  the 
CIF  values  observed  for  fluorescence  time  series  under 

high  hydrodynamic  conditions,  which  cannot  be  distin- 
guished from  the  ClF  values  estimated  for  temperature 

irrespective  of  hydrodynamic  conditions,  and  by  the  ex- 
treme similarity  of  the  empirical  functions  C, T(q ) and 
CF(q)  in  high  hydrodynamic  conditions  (Fig.  7). 

Turbulence  and  biophysical  couplings 

In  order  to  investigate  the  nature  of  the  dependence  be- 
tween temperature  and  phytoplankton  distributions,  we 
subsequently  computed  the  Generalized  Correlation  Func- 
tions, c(p,q)  and  the  related  Generalized  Correlation  Ex- 
ponents, r(p,q),  between  temperature  and  fluorescence 
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time  series,  within  each  subset  defined  above. 
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Figure  8.  The  Generalized  Correlation  Functions  (GCF)  c(p,q ) 
versus  the  time  scale  / in  log-log  plots,  for  (A)  simultaneously 
and  (B)  independently  sampled  temperature  and  in  vivo  fluores- 
cence time  series.  The  function  c(p,q)  shown  here  has  been 
estimated  for  a constant  value  of  the  statistical  order  of  moment 
q of  temperature  fluctuations  (q  = 2),  and  various  values  of  the 
statistical  order  of  moment  p of  in  vivo  fluorescence  (i.e.  p = 
1,  2 and  3,  from  bottom  to  top).  The  slopes  of  the  linear  regres- 
sion estimated  over  the  scaling  ranges  (dashed  lines)  provide 
estimates  of  the  Generalized  Correlation  Exponents  (GCE) 
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Figure  9.  The  Generalized  Correlation  Exponents  (GCE 
r(p,q) , shown  as  a function  of  both  p and  q , which  charac- 
terize in  vivo  fluorescence  and  temperature  fluctuations,  respec- 
tively. The  function  r(p,q)  is  estimated  here  between  time  series 
of  temperature  and  fluorescence  independently  sampled  (C). 

Figure  8 shows  the  GCF,  c(p,q ) , plotted  in  log-log 
plots  versus  the  time  scale  t , for  simultaneously  recorded 
temperature  and  fluorescence  time  series  (Fig.  8a),  as  well 
as  for  temperature  and  fluorescence  time  series  taken  at 
different  moments  of  the  tidal  cycle,  and  a fortiori  inde- 
pendent (Fig.  8b).  Both  the  linear  behavior  of  the  functions 
c(p,q)  over  scales  ranging  between  1 and  20  seconds,  and 
the  positive  values  taken  by  the  GCE,  r(p,q) , estimated  as 
the  slope  of  the  linear  trend  of  the  GCF  c(p,q)  vs.  t in  a 
log-log  plot,  whatever  the  combinations  of  p and  q val- 
ues (Fig.  8a)  confirm  the  validity  of  Eq.  (11).  On  the  other 
hand,  the  weak  values  taken  by  the  functions  c(p,q)  esti- 
mated between  independent  temperature  and  fluorescence 
time  series  (Fig.  8b)  indicate  a low  correlation  between  the 
temperature  and  phytoplankton  biomass  fluctuations, 

((AT,  y ^ and  ((AF,  )p  ^ . This  is  confirmed  by  the  related 

values  of  the  functions  r(p,q) , which  remain  close  to  zero, 
whatever  the  combinations  of  p and  q values  (Fig.  9). 

We  subsequently  refined  these  observations,  comparing 
the  functions  r{p,q)  obtained  between  temperature  and 
fluorescence  time  series  in  different  turbulent  conditions. 
Figure  10  then  shows  the  r(p,q)  obtained  for  all  the  com- 
binations of  p and  q values  with  0.1  increments,  for  three 
levels  of  turbulence  ( 10~4 , 10~5  and  10'6  m2  s'3).  It  appears 
that  the  correlation  between  temperature  and  in  vivo  fluo- 
rescence fluctuations  increases  with  hydrodynamic  condi- 
tions and  is  weaker,  even  nil,  in  low  turbulent  conditions. 

We  confirm  here  the  increased  physical  control  sug- 
gested under  strong  turbulent  conditions  from  the  analysis 
of  the  shape  of  the  function  £,,(<?) , and  its  comparison 
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with  the  function  CrUl)  i see  Figure  7.  On  the  other  hand, 
the  decorrelation  between  temperature  and  phytoplankton 
fluctuations  under  weak  turbulent  conditions  suggests  an 
increase  in  the  biological  contributions  to  the  control  of 
phytoplankton  biomass  distribution,  and  confirms  previous 
observations  (cf.  Fig.  7).  Phytoplankton  distribution  then 
appears  independent  of  the  temperature  distribution  under 
the  lowest  turbulence  levels  investigated  here,  i.e.,  10-4  m2 
s'3  (Fig.  10a). 

More  precisely,  the  shape  of  the  function  r(p,q)  indi- 
cates that  large  phytoplankton  fluctuations  are  associ- 
ated— under  strong  enough  turbulent  conditions — to  weak 
temperature  gradients,  and  conversely.  This  tendency 
seems  to  reflect,  over  a slightly  wider  range  of  scales,  find- 
ings of  Desiderio  et  al.  (1993);  they  observed  the  occur- 
rence of  0.1 -0.2  meter  thick  fluorescent  layers  just  above 
local  temperature  gradients.  While  these  qualitative  com- 
ments suggest  a proportionality  relationship  between  tem- 
perature and  phytoplankton  fluctuations,  the  relation  (16) 
has  never  been  verified.  The  phytoplankton  fluctuations 
then  cannot  be  regarded  as  strictly  proportional  or  ran- 
domly proportional  to  the  fluctuations  of  temperature,  i.e., 
A F = K.AT  and  A F = kAT  , respectively.  In  contrast,  the 
relation  (17)  has  been  verified  testing  the  validity  of  Eq. 
(18)  over  a wide  range  of  b values.  Using  b values  ranging 
between  0.05  and  5 (with  0.05  increments),  we  then 
showed  that  Eq.  (18)  is  verified  for  three  of  the  five  levels 
of  turbulence  investigated  here:  £ = 1 (T4  m2  s"3  with 
6 = 2.69,  £ = 5.1(T4  m2  s'3  with  6 = 2.25,  and  £=10"5  m2 
s'3  with  6 = 1.76.  Figure  11  shows  the  relation  r(p, 0)  vs. 
r(0,bp)  corresponding  to  the  functions  r(p,q)  shown  in 
Fig.  10b  and  c,  with  6 = 2.69  and  6 = 1.76,  respectively. 
The  correlation  shown  between  temperature  and  phyto- 
plankton fluctuations  under  high  turbulent  conditions,  is 
then  related  to  a power-law  dependence  relationship,  i.e. 
A F = k.(AT)1’  . 

Let  us  mention  here  that  at  this  stage  of  the  development 
of  the  multiscaling  techniques  illustrated  in  the  present 
paper,  the  presence  and  the  absence  of  evidence  for  power- 
law  dependence  between  phytoplankton  and  temperature 
fluctuations  cannot  (a  priori)  be  regarded  as  a general  rule 
driving  biophysical  couplings  in  pelagic  ecosystems.  In- 
deed, the  relationship  between  temperature  and  in  vivo 
fluorescence  or,  more  generally,  biophysical  couplings, 
can  be  influenced,  and  then  complexified,  by  both  the 
processes  related  to  phytoplankton  cells  ecology  and  to  the 
sampling  strategy.  Over  the  range  of  scales  investigated  in 
the  present  study,  in  vivo  fluorescence  fluctuations  can 
indeed  result  from  nonlinear  interactions  between  the 
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Figure  10.  The  Generalized  Correlation  Exponents  (GCE 
r{p,q) , shown  as  a function  of  both  p and  q . Here  the  statisti- 
cal orders  of  moment  p and  q characterize  temperature  and  in 
vivo  fluorescence  fluctuations,  respectively.  The  functions 
r(p,q ) correspond  to  three  different  levels  of  turbulence,  i.e.  (A) 

E = IO"6  (B),  10"5 , and  (C)  IO"4  mV3. 


(shown  recently  to  influence  both  spectral  and  multifractal 
properties  of  phytoplankton  biomass;  Seuront,  1999; 
Lovejoy  et  al.,  2001),  and  the  turbulent  processes.  The  re- 


physiological  state  of  phytoplankton  cells,  their  related  suiting  fluctuations  can  then  be  potentially  more  complex 

coagulation  ability,  the  specific  composition  of  phyto-  than  those  observed  for  a purely  passive  scalar  such  as 


plankton  assemblages,  the  grazing  impact  of  copepods  temperature. 
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Figure  11.  Plots  of  the  Generalized  Correlation  Exponents 
r{p, 0)  versus  r(0.bp)  for  two  levels  of  turbulence,  £ = 1 0~ 

with  6 = 2.69  (A)  and  £ = 1CT'  with  6 = 1.76  (B),  The  slopes 
of  the  regression  lines  (continuous  line)  cannot  be  distinguished 
from  that  of  the  relation  bp  = q (dashed  line).  This  shows  the 
validity  of  Eq.  (18)  and  the  symmetry  of  the  functions  r(p,q)  in 
the  bp-q  plane. 

On  the  other  hand,  considering  (i)  the  highly  intrusive 
character  of  the  sampling  device  used  to  record  tempera- 
ture and  in  vivo  fluorescence  (a  Sea-Bird  25  Sealogger 
CTD  probe  and  a Sea  Tech  fluorometer,  respectively),  and 
(ii)  the  spatial  separation  of  the  temperature  and  fluores- 
cence sensors  (=10  cm),  it  is  reasonable  to  think  that  the 
most  extreme  events,  which  are  also  the  most  infrequent 
and  contribute  to  the  highest  values  of  the  moments  p and 


q , are  not  necessarily  simultaneously  sampled  in  a proper 
way,  leading  to  spurious  joint  high-order  properties.  In 
particular,  this  could  explain  the  weakening  of  the  Gener- 
alized Correlation  Exponents,  r{p,q),  for  values  of  the 
moments  p and  q such  as  p > 4 and  q > 4 . The  latter 
limitation  will  be  easily  fixed  in  the  near  future  with  the 
increasing  development  and  availability  of  microstructure 
profilers  able  to  record  both  physical  and  biological  pa- 
rameters with  a high-spatial  resolution  from  mildly  intru- 
sive sensors  ( Seuront  el  a /.,  200 Id, e).  The  former  limita- 
tion will  be  difficult  to  overcome,  because  of  the  intrinsic 
physiological,  biological,  and  ecological  processes  affect- 
ing phytoplankton  populations. 

Conclusions:  the  cost  and  gain  of  increasing 
analysis  complexity 

Our  results  demonstrate  that  single  scaling  properties  of 
temperature  and  phytoplankton  would  have  erroneously 
led  to  conclude  that  the  latter  could  be  regarded  as  a purely 
passive  scalar.  However,  investigations  of  their  detailed 
variability  in  the  universal  multifractal  framework  lead  to 
consideration  that  phytoplankton  biomass  presents  very 
specific  properties,  depending  on  the  intensity  of  local 
turbulent  processes.  This  suggests  a differential  control  of 
phytoplankton  biomass  distribution  involving  complex 
interactions  between  turbulent  velocity  fields  and  the  in- 
volved particles  properties  such  as  buoyancy,  size,  density 
or  aggregative  properties.  In  particular,  this  assertion  is 
supported  by  the  specific  behavior  of  phytoplankton  mul- 
tiscaling properties  under  low  turbulent  conditions,  sug- 
gesting the  prevalence  of  phytoplankton  cells  specific 
properties  on  turbulent  processes. 

This  has  been  confirmed  and  specified  by  the  introduc- 
tion of  the  Generalized  Correlation  Functions/Exponents. 
The  correlation— and  then  the  couplings — between  tem- 
perature and  phytoplankton  biomass  fluctuations  are  thus 
higher  under  high  hydrodynamic  forcing  and  decrease  pro- 
gressively to  collapse  under  low  hydrodynamic  forcing.  In 
particular,  careful  examination  of  the  functions  r{p,q) 
indicates  that  large  phytoplankton  fluctuations  are  associ- 
ated— under  strong  enough  turbulent  conditions — with 
weak  temperature  gradients,  and  conversely.  These  results 
are  fully  congruent  with  more  intuitive  and  qualitative  re- 
sults by  Desiderio  et  al.  (1993)  over  similar  scales  and 
Wolk  et  al.  (2001)  at  smaller  scales.  Whatever  that  may  be, 
even  if  additional  work  and  data  sets  will  be  needed  to 
achieve  more  definitive  functional  scenarios,  the  joint  mul- 
tifractal framework  introduced  here  provides  the  first  ob- 
jective method  to  characterize,  both  qualitatively  and 
quantitatively,  the  details  of  the  relationships  between  any 
pair  of  simultaneously  recorded  patterns  or  processes.  In 
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particular,  this  framework  makes  no  assumptions  about  the 
spectrum  or  the  distribution  of  either  data  sets,  and  it  takes 
fully  into  account  their  intrinsic  multiscaling  properties. 
We  should  also  note  here  that  this  work  provides  the  very 
first  evidence  of  the  multifractal  character  of  ocean  turbu- 
lent velocity  fluctuations,  and  reveals  their  extreme  simi- 
larity in  the  ocean  and  the  atmosphere. 

Finally,  regarding  the  additional  quantity  of  information 
brought  by  the  univariate  and  bivariate  multiscaling  proce- 
dures described  here  to  our  understanding  of  the  structure 
of  velocity,  temperature  and  phytoplankton  fluctuations,  as 
well  as  of  the  couplings  between  temperature  and  phyto- 
plankton biomass,  we  claimed,  as  already  stressed  else- 
where (e.g.,  Yamazaki,  1993;  Seuront,  2001 ; Seuront  et  al., 
2001c),  that  taking  into  account  the  reality  of  the  intermit- 
tent nature  of  bio-physical  microstructures  could  be  the 
first  step  towards  a general  understanding  of  structures  and 
functions  in  marine  pelagic  ecosystems.  The  quite  elevated 
cost  associated  with  making  use  of  the  single  and  joint 
multifractal  tools  presented  here  in  terms  of  algorithmic 
developments  and  computation  time  consumption  will  then 
be  easily  overcome  by  the  priceless  gain  in  emerging 
knowledge. 
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